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q . Abstract 

Nuclear anapole moments of 133 Cs, 203 > 205 T1, 207 Pb, 209 Bi are treated in the 
single-particle approximation. Analytical results are obtained for the oscillator potential 
without spin-orbit interaction. Then the anapole moments are calculated numerically 
in a Woods-Saxon potential which includes spin-orbit interaction. The results obtained 
demonstrate a remarkable stability of nuclear anapole moment calculations in the single- 
particle approximation. 
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1 Introduction 

The existence of parity nonconservation (PNC) in atoms is firmly established at present (see, 
e.g., Ref. [[T|). To be precise, only the nuclear-spin- independent PNC effects in heavy atoms 
have been observed up to now. Just these effects are enhanced as Z 2 Q. The last enhancement 
factor, the so-called weak nuclear charge Q which is numerically close to the neutron number 
N = A — Z, is due to the fact that in the nuclear- spin- independent phenomena all the 
nucleons act coherently. 

As to the atomic PNC effects dependent on nuclear spin, they evidently lack this coherent 
enhancement and are therefore much smaller. There are strong reasons to expect that these 
effects are dominated by contact electromagnetic interaction of electrons with nuclear anapole 
moment (AM) g §. 

Anapole is a new electromagnetic moment arising in a system without centre of inversion 
M. It exists even in such a common object as a chiral molecule in a state with nonvanishing 
angular momentum ||. Nuclear anapole moment is induced by PNC nuclear forces. 

The electromagnetic PNC interaction of electrons with nuclear AM is conveniently char- 
acterized in the units of the Fermi weak interaction constant G = 1.027 x 10~ 5 m -2 [] by a 
dimensionless parameter k, (Its definition is given in the next section.) A closed analyti- 
cal expression has been obtained for this constant || within the nuclear shell model under 
some extra simplifying assumptions. In particular, for 133 Cs this model prediction is (at the 
numerical value of the effective PNC nuclear constant accepted in the paper quoted) 

K ( 133 Cs) = 0.33. (1) 

More reliable numerical calculations with the Woods-Saxon potential including the spin-orbit 
interaction result at the same PNC nuclear constant in |3j 

K ( 133 Cs) = 0.25. (2) 

The constant k( 133 Cs) has been calculated also in Ref. . In that paper only the 7r-meson- 
exchange contribution to the P-odd nuclear forces has been included. This contribution 
constitutes roughly a half of the strength of these forces according to the estimates accepted 
in Ref. || (the second half is due mainly to the P-odd p-exchange). No wonder therefore that 
the number 0.14 obtained in Ref. || for k( 133 Cs) constitutes roughly a half of the total result 
(0). (In Table [l], containing the summary of theoretical results, in the corresponding entry 
we indicate in brackets what to our guess would be the result of Ref. || if the p-exchange 
were included.) 

Later numerical calculations with the oscillator potential |7|] led at the same value of the 
P-odd nuclear constant g p to 

K ( 133 Cs) = 0.24. (3) 

Less satisfactory is the agreement between the results for the AM of 209 Bi obtained in Refs. 
0, 0. As to the AM of 203 > 205 T1 the disagreement between the predictions made in Refs. 
||. H is even stronger. (The values of k( 203,205 T1), K( 209 Bi) at the same magnitude of g p as 

4 The system of units with % = 1, c = 1 is used; m is the proton mass. 



accepted in Ref. || is not explicitly given in Refs. []7|, ||. So, in the corresponding entries of 
Table [I] we present in brackets our extrapolation of k( 203 ' 205 T1), K( 209 Bi), obtained in Refs. 
|7|. Pj, from their values of g p to that accepted in Ref. 0.) 

Experiments aimed at the detection of nuclear AM in cesium, thallium, lead, bismuth are 
underway in many groups. The first evidence of a nuclear-spin-dependent P-odd effect has 
been seen already in cesium 0. The result of this experiment is 

K ( 133 Cs) = 0.72(39). (4) 

Therefore, detailed theoretical investigation of various contributions to nuclear AM looks 
quite relevant. 

The present paper is organized as follows. In the next section we describe the constant- 
core-density approximation which leads to a closed analytical expression for nuclear AM, 
and introduce some necessary notions and relations. Then we consider somewhat more 
sophisticated model which still allows for an analytical treatment. This model gives some 
idea about the corrections to the above leading approximation. But the most essential part of 
the paper is devoted to serious numerical calculations using a realistic description of the core 
density. The nucleon wave functions and Green's functions are obtained with a Woods-Saxon 
potential which includes the spin-orbit interaction. The contribution of the current generated 
by the spin-orbit interaction (omitted in Ref. 0) is taken into account. As distinct from 
Refs. [0, Q, the spin-orbit interaction itself is treated beyond the perturbation theory. 

Though the deviations of the contact and spin-orbit currents from naive potential ex- 
pressions are included, but in all other respects we restrict throughout the present paper to 
the single-particle approximation, that of a valence nucleon above a spherically-symmetrical 
core. The many-body effects are certainly of importance for nuclear anapole moments. We 
shall consider them in future publications. 

2 The leading approximation 
for nuclear anapole moment 

PNC interaction in a system (atomic nucleus is now of interest to us) mixes opposite parity 
states of the same total angular momentum and creates in it a spin helical structure |l], @j]. In 
this way such a system with nonvanishing magnetic moment acquires a specific configuration 
of magnetic field, of the type created by the toroidal winding. This is what is called anapole 

E- 

The AM vector can be conveniently defined as [[l], [| [| 

a = —n / dfr 2 j{f) (5) 



where j(r) is the current density operator. The vector-potential produced by the AM is 

A(f) = -[oA-V(oV)]- >aS(f). (6) 



We omit in the last expression the term V(a V)l/(47rr), which can be obviously eliminated 
by a gauge transformation. 

When calculating the AM of a heavy nucleus we restrict ourselves to the shell model and 
one-particle approximation (both in what we will call the leading approximation and beyond 
it). We shall assume that the nuclear spin / coincides with the total angular momentum of an 
odd valence nucleon, while the other nucleons form a core with the zero angular momentum. 
The effective P-odd potential for an external nucleon can be presented as follows: 

w = 71 2m~ d ^p^ + p( r )p]- ( 7 ) 

Here a and p are respectively spin and momentum operators of the valence nucleon, p(r) is 
the density of nucleons in the core normalized by the condition / dfp(r) = A (the atomic 
number is assumed to be large, A ^> 1). The numerical value of the dimensionless constant 
g p in the case of an external proton is perhaps close to 4 — 5 (see. below). For an external 
neutron the corresponding constant g n is smaller, most probably g n <l. 

The leading approximation for the AM of a heavy nucleus corresponds to the assumption 
that the density p(r) is constant throughout the space and coincides with the mean nuclear 



density po- This approximation, first used in Ref. flQ |, is reasonable if the wave function of 
the external nucleon is mainly localized in the region of the core. The Schrodinger equation 
for the external nucleon 

[ — —A + U(r) + W(f)Wf) = Eij)(i f ) (8) 

2m 

to first order in W for p(r) = Po = const has the elementary solution 

G 
if){f) = (1 - i-^gp ar)ip {r)- (9) 

Here ipo(r ) is the unperturbed wave function of the external nucleon. It might seem that now 
interaction ([7|), which is equivalent to the electromagnetic interaction with a constant vector- 
potential A = —(G/y/2)(g/e)p a, should not result in any physical effects at all. However, 
the spin part of the current density 

Hn = |^-tf x ww,), (io) 

2m 

(/i is the nucleon magnetic moment) does work due to the noncommutativity of the a- 
matrices, even in this approximation. Simple calculations using formulas (|5|), @ and ( ]T0|) 
yield 



Ggpo 2nep . 2 . KI 



V2 m N '/(/ + 1)' 

K = (I + l/2)(-) I+1/2 - 1 . (11) 

Here I is the orbital angular momentum of the external nucleon. As to its mean square radius 
(r 2 ), it coincides to good accuracy with the squared charge radius of the nucleus 

r 2 = (3/5)R 2 = (3/5)r 2 A 2 / 3 , r = 1.2 fm. (12) 



It is useful to present also the effective local AM operator which, acting in the space of 
nonperturbed wave functions ipo(r), produces the result ([11]). This operator is 

- Ggpo 2vre/i _> 2 _,^ni n „\ 

-or + r [or )\. (13) 



\/2 m 
Setting p = (47rro/3) _1 , we finally obtain from (|TT 



G 9 en .o/q KI 

a = —=—g^—A 2 ' 3 — -. (14) 

y/2l0mr 1(1+1) K ' 

The A-dependence of the AM is very natural. Indeed, since the anapole corresponds to 
the magnetic field configuration induced by a toroidal winding, the AM value should be 
proportional to the magnetic flux, i.e., to the cross-section area of the torus. This is the 
origin of (r 2 ) in formula (|TTD and of A 2 !' 3 in (|j~4]). 

Let us turn now to the PNC problem in atoms. The Hamiltonian of the interaction of an 
electron with vector-potential (^) can be presented as 

H a = eaa6(f ) = — K S(f ) (15) 

(— e is the electron charge, a are the Dirac matrices). The Fermi constant G serves as 
the natural unit for the AM, that arises in first order in the weak interaction, and has the 
dimension cm 2 . In this unit a convenient characteristic of the nuclear AM for the atomic 
PNC problem is the dimensionless constant k. According to (TTJ]) it equals || 

k=^— A2 "- ( 16 ) 

10 mro 

The enhancement ~ A 2 ' 3 compensates to a large degree the small fine structure constant 
a = 1/137. That is why the nuclear AM is perhaps the main source of the nuclear-spin- 
dependent PNC effects in heavy atoms @, |3|. The constants k for nuclei of experimental 
interest, as given by simple analytical formula (|i~6|) , are presented in Table |l[ Their numerical 
values for Cs, Tl, Bi correspond to g p = 4, the number assumed in Ref. 0. In the next 
line of Table [l] the results of numerical calculations || are presented at the same g p . Those 
calculations were performed using a realistic description of the core density p(r). The wave 
function and the Green's function of the valence nucleon were calculated with a Woods- 
Saxon potential which included the spin-orbit interaction. The crude analytical calculation 
is obviously in reasonable agreement with that numerical one. 

3 Analytical treatment of nuclear AM in a single par- 
ticle oscillator potential. 
Contact current contribution 

Our analysis of various contributions to nuclear AM beyond the leading approximation we 
will start from a more realistic model which still allows for an analytical treatment. The 



model consists in the use of the oscillator potential for the valence nucleon and in neglect of 
the spin-orbit interaction. 

Let us consider at first the spin current contribution to the AM. Substituting expression 
(p"0|) into formula (|5|) and integrating by parts, we transform the corresponding AM operator 
to the following form: 

a s = — -fxa. (17) 

m 

When using the oscillator potential 



U{r) 






2 
for the valence nucleon, its radius-vector transforms to 



mu) 2 



[H,p\ (19) 



where H is the valence nucleon Hamiltonian. Substituting this expression into the standard 
second-order perturbation formula for AM 

vre/i ^ (0\[H,p x a]\n)(n\W\0) + (0\W\n)(n\[H,p x a] 



-i- 



\p V J \\. J - L il J ^ "Jl"'/\"'l vv l w / ~ \ yJ \ VY \' b l\"'\[- L± iF 2 U W^I /OQ^ 



m 2 uj 2 V E - E n 

we reduce it by means of the completeness relation to 

a s = -z^(0|[px«7,W]|0> (21) 

With the explicit form (0) for the weak interaction Hamiltonian W this formula can be 
rewritten after elementary transformations as 

a, = % ^TT^TT 4^<0|PP 2 - W - '&m- (22) 

V2 m 1(1 + 1) m l uj A 2 2Kr 



Let us note here that previously the oscillator-potential model was used in Ref. 11 for 



the investigation of other P-odd nuclear effects. Then it was applied in Ref. [[12] to estimate 



the contribution of core excitations and to rederive formula ( ^4|) (which follows indeed from 
expression (|22D at p(r) = po = const when taking into account the oscillator-potential version 
of the virial theorem: (0|p 2 /m|0) = (0\mu 2 r 2 \0)). 

We pass over now to the orbital contribution to the AM. As it was demonstrated in Ref. 
0], the AM of a charged particle in a spherically symmetrical potential can pe presented in 
the following form: 

a = -^(/i - |)j7J^i) E Won - y (0\r 2 j° c - r(rj° c )\0). (23) 

Here r]Q n is the P-odd admixture of an intermediate state \n) to the initial one |0) (this 
coefficient is purely imaginary), ro n is the matrix element of r between those states, j® is the 
contact current operator: 

j° c =te[W,f] = ^-pa. (24) 

V2m 



Therefore, to take into account within the simple-minded potential approach the orbital 
contribution for a proton we should substitute /i — 1/3 for /x in formula fl22|) and add to that 
expression the following contact term: 

Gg2ixe KI 1 /nl 2 



^2 m 1(1+1)3 



<0|rV|0). (25) 



In the case of a valence neutron formula ( p2]) needs obviously no modification within the 
same model. 

However, in fact the contact current contribution both for a valence proton and neutron 
looks differently. Indeed, let us start from the weak interaction Hamiltonian for the outer 
nucleon a(a = p,n) constructed from corresponding two-body operators (see, e.g., book [[J): 

C 1 

W {2) = —S J— £ ({(Sab&a - 9baVb) ' (fa - Pb) , 5{r a - f%) } + g' ah \o a X S h ] • V<J(f a - ffc)) , 

V ^ a,fe 

(26) 
where the notation { , } means anticommutator. After averaging this expression over the 
core nucleons we obtain formula (0) with 

Z N 

9 = 9ap^+9an^- (27) 



The "best values" [13] of weak coupling constants in the one-meson-exchange approximation 



lead to the following numbers for the constants discussed H [14], |l5f| : 

9pp = 9nn — t.o, g pn = b.o, g n p = Z.Z, \^°) 

and correspondingly to 

9 = 9 V = 4-5 (29) 

for the valence proton in Cs, Tl, Bi, and 

9 = 9n < 1 (30) 

for the valence neutron in 207 Pb. Here the constants g pp , g nn , g p and g n are effective ones, 
they include already the exchange terms for identical nucleons. These constants include also 
additional suppression factors reflecting long-range and exchange nature of the P-odd one- 
meson exchange as well as the short-range nucleon-nucleon repulsion. Their values agree with 
those obtained in Refs. [|17], |18| and differ from those given in Ref. |19] since we include into 



the p-meson exchange the suppression due to the short-range repulsion. 

The electromagnetic interaction can be introduced by changing p =>■ p — eA(r). The 
current density operator is then a derivative of a Hamiltonian over A(f). This is equivalent 
to change one power of a momentum operator p a =^ e a 5(r — r a ) in Hamiltonian. This 
procedure leads to the following expression for the contact current density operator: 

I = ~{£[^Var„U(r-r a )} = E^W 7 " ^), (31) 



2 



where 

- G 1 — 

Jc(r) = -7^ — &a X]( e « _ Cb)gabS(r- T a ). (32) 

For the valence protons e a = e and this operator reduces to 

Jc = —7=—Og pn J2 S (^ ~ f «)- ( 33 ) 

For the outer neutron (e a = 0) we get 

^-Jl^P^-^' (34) 

After averaging formulae (^), (0) over the core nucleons we get the following effective 
contact currents for the valence proton and neutron: 

Jc = ~ftm 9pn ^ P ^ = ^ + ^' ^ 

G e Z 

3c = -7i™ 9pp A pir)a; 

1: = -^9n P j P (r)a. (36) 

They differ obviously from the naive ones, j° for a valence proton (see (|24D ) and zero for 
a valence neutron. Technically it is convenient to retain previous, naive results for anapole 
moments, supplementing them with the following correction terms: 

3* = -7r(0|r 2 S|0>, (37) 

S" 1 = -^(0|r 2 J:|0) (38) 

for valence proton and neutron respectively. In this way we get the following closed expression 
for the AM of a nucleus with a valence proton: 

G 2ne KI , 2v 



v/2 m I(/ + l) 
■ (0|^[2p-rV(O- 2 _ ( ; +3/2) pX^ + ^)] 

->- 2 ^> - 2^+3/2) ^ + ^ )]} (39) 

- ^(1-^)^10); 



<£ = g pp Z/A = 0.6, (r 2 ) = (N + 3/2)/mu> 
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where N is the oscillator principal quantum number. In the case of a valence neutron the 
result is 

G 2vre KI . 2v 



y/2 m 1(1+ 1) 

■ (0\g n fi n [2p - r 2 p/(r 2 ) —^ — p'(d r + l !±^- )] 

\ \yn^ni h hi\ I 2 mu{N + 3/2) H v Kr n 

- ^n(l-^)r 2 p|0); (40) 

g 1 n = g np N/A = -1.3. 

To get the numerical values of nuclear AMs we assume for the core density a step-like 
profile f(r): 

p(r)=p f(r) = p 9(R-r). (41) 

As to the mean square radius of the valence nucleon (r 2 ), it is natural to identify it in 
the shell model with the nuclear magnetic mean square radius (r^). The empirical data 
on the latter (referring unfortunately to lighter nuclei only) are presented in Table |2|. The 
observation is that in the 41 Ca2o, with a valence neutron, (r, 2 ^) is close within the error 
bars to the value predicted by the oscillator model. Meanwhile for nuclei with odd proton, 
45 Sc2i, 51 V23, 59 Co27, its value is much better approximated by the naive formula (]I2"D. 

For further calculations it is natural to go over to the usual dimensionless oscillator 
variable PU 

x = r 2 lp\ p 2 = (me)" 1 = -f-r'rtA^. (42) 

5 3 

Then the expectation values entering eqs. (|3~9"D, ( fiOD reduce to 

(0|/(r)|0) = T ^ { t+!/2) £ d ^- X * l+1/2F2 (- n > l + 3 A x), (43) 

(0|r 2 /(r) |0) = P ^t^ljl] J^ dxe-*x l ^F*(-n, I + 3/2, x), (44) 

(o\ f'( r )(d + l ( l + 1 h \0) = 2 r(n + Z + 3/2) x/2 {l+3)/2 

2mw(N + 3/2) {(,lJ[ )[ ° r+ Kr m (N + 3/2)n!r 2 (Z + 3/2) C 

■F(-n, I + 3/2, X)( ± + l l±^) e -^X^F(-n, I + 3/2, X). (45) 

Here X = R 2 /p 2 = (5/4)(3/2) 1/3 A 1/3 = 1.43A 1/3 , F(a,b,x) is a degenerate hypergeometric 
function, n= (N -l)/2. 

We are interested in the following nuclei: 



Cs 1 = 7/2, 1 = 4, K = 4, N = 4, n = 0, X = 7.304; 

Tl I = 1/2, 1 = 0, K = -1, X = 4, n = 2, X = 8.423; 

Bi J = 9/2, / = 5, X = 5, N = 5, n = 0, X = 8.492; 

207 Pb J = 1/2, / = 1, if = 1, N = 5, n = 2, X = 8.464; 
For Cs and Bi where n = the integrals entering expectation values fl4"3]), (f4"4|) are 



nothing else but well-known (and tabulated, see, e.g., pl| ) incomplete gamma-functions, 
T(7 + 3/2, X), r(Z + 5/2, X) respectively. In the case of Tl and Pb those expectation values 
can be also reduced to incomplete gamma-functions. 

The numerical results for nuclear AMs, obtained in this model, are presented in Table |l|. 
For Cs, Tl and Bi, at the same value of g p , they are about 10% smaller than those of the 
leading approximation. The correction to the spin current contribution constitutes about 
-4% in Cs and Bi, -7% in Tl. The total contribution of the convection and contact currents 
is negative in cesium, thallium and bismuth. As to lead, the contact current contribution, 
though nonvanishing, is very small numerically. 

4 Nuclear AM in Woods-Saxon potential. 
Spin-orbit current contribution 

Though being convenient for analytical treatment, the step-like density profile and oscillator 
model for a single particle potential are however too crude. Moreover, the above analytical 
model does not take into account the spin-orbit interaction which is, as we will see below, quite 
essential for numerical results. So, in this section we will describe a numerical treatment of a 
much more realistic description of a nucleus based on the Woods-Saxon potential including 
spin-orbit interaction and on a more realistic description of nuclear density. 

The profiles of both density and the central part of nuclear potential are known to be 
similar and well described by a Fermi-type function 

1 + exp( 1 ^) ' 
So, the total single-particle potential U(r) has been chosen in a standard Woods-Saxon form 

U(r) = U f(r) + U ls -^P-{W) + U c {r), (47) 

r dr 

where Uc{r) is the Coulomb potential of a uniformly charged sphere. In order to study 
stability of AM calculations against variations of the single particle potentials we shall use 
several sets of the Woods-Saxon potential parameters 

9 



The solution of eq. (§) can presented in a form similar to (|9|) 

W) = ^o(r) + <ty(f) = ^o(r) - t^=g Po (an) SR(r)) n Ilm (n), (48) 

where Qu m (n) is a spherical spinor, n = f/r. The correction Sip is of the parity opposite 
to that of the initial state ipo. The radial function SR(r) is normalized in such a way that 
without spin-orbit potential and for constant density it is (see (||])) 

SR(r)=rR (r). (49) 

This correction can be expressed via the Green's function Gn'{r, r') of the unperturbed radial 
Schrodinger equation for the orbital angular momentum V = 21 — I: 

SR(r) =~J r ,2 dr' [±G IV {ry)R {r') - Gll ,(r,r')-^R (r') 

*G IV {ry)R Q {r')}f{r>). (50) 



IK 



After presenting the Green's function through two linearly independent solutions of the radial 
Schrodinger equation we get 

5R(r) = ^1 f dr' (u$(0«o(0 - uf> '(r'H(r') + ^<4V)«o(r')) f(r') 

, U fk r ) r , / / (1)/ A //A (1)// A , A , 2 ^ 



r Jo 



dr f uSWoO"') - u^'(r')uo{r') + —u$ (r')uo(r')) f(r'), (51) 



where Uo(r) = rRo(r); u^(r) and u^(r) are solutions regular at the origin and at the infinity 
respectively. Those last two solutions are normalized to the unit Wronskian: vS^u^' — 
uM'u® = 1. 

The results of this calculation of the 5R for lh 9 / 2 proton state in 209 Bi together with the 



leading approximation ( ^9|) are shown in Fig.l. In the same picture we demonstrate how SR 
is influenced by the spin-orbit potential, as well as by the deviation of the density from a 
constant one. The real form of 5R differs considerably from the leading approximation, due 
mainly to the spin-orbit interaction. As to the realistic density profile, it makes SR slightly 
smaller as compared to the constant density approximation. 

All contributions to the anapole moment can be expressed via the radial correction SR. 
The spin current term can be obtained from eq. (p|) by substituting 

/•CO 

{SR\r\R Q ) = / r 2 dr SR{r)rR Q {r). (52) 

Jo 

for (r 2 ). So, this contribution to the dimensionless anapole constant equals 

k s = 27rg^^(SR\r\R ). (53) 

m 

10 



To obtain the convection current contribution we substitute corresponding current density 
operator 



into formula (|5|) which gives 



-^£{V^(r--f p )} (54) 



-rr^KAr 2 }!^) (55) 

m 



and 

n conv = -7rg^(5R\r 2 (^- + ^^)\R ). (56) 

mK dr r 

When calculating the contact current contributions it is convenient again to split it into 
the "naive" part arising from a single particle potential (and vanishing for a valence neutron) 
and the part which is due to the current densities fl35"D, fl36|). 

The contribution of the total contact current (|35D to the AM can be expressed via radial 
matrix element as 

k c = k° c + kI < = 71^(1 - ^)(r 2 /(r)), i = 0, 1 (57) 

m 2K 

where g® = g p , g® = 0, and g* were defined in fl39l), (|40|). 

One more term in AM originates from the momentum dependence of spin-orbit interac- 
tion. As it will be seen below, this is the most significant correction to the leading approx- 
imation. The spin-orbit term in a single particle potential (f|7|) for a proton generates the 
electromagnetic current density 

jf s = ie[U(f),r\ = eU ls ^-±a x n. (58) 

For an outer neutron the corresponding current density vanishes. The spin-orbit contribution 
to the proton AM and to k is respectively: 

4 = 27T-^=egp U ls (8R\r 2 f(n(an) - a)\R ) = 2n-^=egp U ls (SR\r 2 f\R ) ™ , (59) 

K ls = 2nagp Q U ls (5R\r 2 f\R ). (60) 

However, as it was the case with the contact current, this is not the complete result. The 
true spin-orbit current must be obtained from two-particle spin-orbit interaction which can 
be written as 

Uif = \Y, < (Pa - p b ) ■ (*« + a b ) x VS(r a - r b ). (61) 

ab 

Averaging it over nuclear core we obtain the spin-orbit part of the single-particle potential 
(|7|) with 

Uis = (itfj + UlTj) Po- (62) 

11 



The spin-orbit interaction constants were fitted in Ref. p6| : 



U[ s p = U™ = 36.6 MeV ■ fm 5 , Uf s n = U%> = 134.3 MeV ■ fm 5 . (63) 



A word of caution is proper here. We will understand expression (|61|) as a local limit for an 
operator of a nucleon-nucleon interaction with the proper tensor structure. For the identical 
nucleons the exchange interaction doubles the direct one. The values of Uf s p , Up s n presented 
in (|63| ) are just those doubled, effective constants. 

The spin-orbit current density from the interaction (|6l|) is 



3u = i E Ms\ Za?a] S(f- r a ) = £ ug e a 5(f - f a ) (a a + a b ) x V5(f- f b ). (64) 

a ab 

This expression for the two-particle spin-orbit current density was used earlier in the de- 



scription of magnetic properties of a nucleus in ||24|| . Including the spin-orbit current into 
the theoretical description of nuclear magnetic moments and magnetic transition amplitudes 
improves the agreement with experiment. 

In the current density generated by the pp part of operator (|6T| ) the exchange and direct 
contributions cancel each other, i.e., the contact spin-orbit pp interaction does not generate 
a current at all. Then, averaging expression (|64]) over the core nucleons and separating the 
single particle contribution (|58|) , as in the case of the contact current, we get the following 
effective spin-orbit current for the valence proton: 

> TTP n ™ df(r) _, ->q -^ 

Jis = eUf s Pv-^—^r<? x n = j l8 + j ls ; (65) 

^i ttvv Zdf(r)_ 

] >- = -*""°i * " x "■ 

The correction j} s to the "naive" spin-orbit proton current is relatively small being suppressed 
by a factor 

K ~ zuf! + NUf s n 

as compared to the potential part. This factor constitutes -0.16 for cesium and -0.15 for 
thallium and bismuth. 

It is noteworthy that even for a valence neutron the spin-orbit current does not vanish. 
It is generated by the neutron-proton interaction and equals 

lis = ~eU??j p(r) V x (r(f)a^(r)) (66) 

This expression is similar to the spin current density ( |i~0l) and renormalizes in an obvious 
way the magnetic moment of an outer neutron. 

It is convenient to single out in the sum of all contributions to the effective constant k 

K = K s + Ki s t Ki s ~r Kconv > ^c ' ^c = ^ ' ^ls ' ^ci \® ' ) 
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the "naive" potential one k°. The corrections to it k} s and k\ in fact go beyond the single- 
particle approximation. The single-particle contribution can be written in a more compact 
form, analogous to (E3|): 

3° = (1 - ^K + \< - y (0|r 2 £° - r(f£o)|0). (68) 

This expression was used to check the accuracy of our numerical calculations which turned 
out very high. 

To study the stability of AM calculations under the variations of the single-particle poten- 
tial we performed them for two sets of the Woods-Saxon potential parameters (f|7]). In both 
sets the radii and diffuseness parameters, R and a, are the same for the profiles of central 
and spin-orbit parts of potential. 

The first set is p0|: 



R = R ls = 1.27A 1/3 fm, a = a ls = 0.67 fm, 

N - Z 
U = (-51 ± 33 — - — ) MeV, U ls = -0.35 U Q , (69) 

where the signs + and — refer to neutrons and protons respectively. 

The second set p5j has somewhat different parameters, but gives as good fit to the single- 



particle level positions. It is: 



R = R ls = 1.2AA 1/3 fm, a = a ls = 0.63 fm, 



N - Z N - Z 

C/o = (-53.3 ± 33.6 — — ) MeV, U ls = -0.263 (1 + 2 — j—)U (70) 



In both cases we assume the same values of the density parameters pO 



R = 1.11 A l/3 fm; a = 0.54 fm p = 0.17 fm~ 3 . (71) 

The results of calculations obtained for both sets are listed in Table |^. The main contri- 
bution to AM, as expected, comes from the spin current. The variation of the parameters of 
potential changes this term by 3 - 6%. The next in magnitude term is due to the potential 
part of spin-orbit current, it constitutes about 30% of the spin current contribution. 

In the calculations based on the sets of parameters (|69|) and (|70|) the LS-potential was 
not consistent with the two-particle LS-interaction (|61|). So, the potential LS-contribution to 
the AM was calculated in two ways. First, starting from the LS-potential itself. Second, it 
was obtained from the spin-orbit potential constructed from the two-particle interaction (|6T| ) 
- (^3]). The difference between second and first ways of calculation is listed as Sk^ in Table |3|. 
It is only natural to perform a consistent calculation, choosing the amplitude of LS-potential 
in accordance with (0), (|63|) , and its radius and diffuseness parameters the same as those of 
nuclear density. In this way we come to the third, LS-consistent, version of the potential: 

TV — Z 
R = 1.21 A 113 fm, a = 0.67 fm, U = (-51 ± 33 — - — ) MeV, 

N — Z 

R ls = l.llA 1 / 3 fm, a ls = 0Mfm, U ls = (14.5 =f 8.3 ~^~) MeV. (72) 
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We kept the central potential parameters as in (|69| ) because it gave smaller Sk^ s . The calcu- 
lation based on this version of the potential give Sk^ s = 0. 

It is noteworthy that the sum of all contributions to AM is less sensitive to a specific 
potential than each of them taken separately. Its variation does not exceed 5% when going 
from one set of parameters to another. 

We believe that the results obtained with the LS-consistent potential are the most reliable 
ones. They are presented in the last line of Table |l|. 

5 Discussion of the results 

Let us discuss now why the present results differ from those obtained in Refs. |3], [7|, || . 

As distinct from Ref . || , we have included now the contribution of the spin-orbit current 
which turns out quite essential. 

On the other hand, the Woods-Saxon potential used here for numerical calculations is 
more realistic than the oscillator one used for those calculations in Refs. 0, ||. Then, in 
Refs. |7|, H the spin-orbit interaction is treated perturbatively, while in the present work 
it is treated exactly. Let us note here again that the oscillator potential without spin-orbit 
interaction allows one to get an exact analytic solution for the AM, as it has been done in 
Section 3 of the present work. 

To summarize the comparison we wish to say that the above arguments can serve as 
serious grounds to consider the numerical results of the present work to be the most reliable 
ones for the single-particle approximation. 

However, accurate quantitative predictions for nuclear AM's cannot be made without 
proper treatment of nuclear many-body effects. As to the core excitations by the weak 
interaction Hamiltonian (pS|) , their contribution to nuclear AM has been demonstrated to 
be small [0. However, the problem of the configurations mixing caused by usual P-even 
residual nucleon-nucleon interaction is here much more serious than in the case of nuclear 
magnetic moments (or the second neutral current constant 0). The point is that in the 
last cases it is the specific, kinematical, nature of the operator a which allows one to restrict 
with a reasonable accuracy to the excitations within a spin-orbit doublet only. As to the 
AM, even if the multiconfiguration problem could be reduced to the calculations with the 
effective operator fll3|) (which demands a special proof by itself), that operator is coordinate- 
dependent and in this sense resembles more the magnetic octupole operator than the magnetic 
dipole one. And many-body effects can renormalize M3 operator much more strongly than 
Ml operator. A drastic example in this respect is 209 Bi where the one-particle value of the 



magnetic octupole moment is almost four times smaller than the experimental one |24]. 

Let us note here that due to the theoretical uncertainties discussed, in the situation when 
the AM contribution dominates essentially atomic nuclear-spin-dependent PNC effects, the 
proposal |7| to single out in those effects the true neutral current contribution by combining 
experimental data from various heavy nuclei, does not look realistic. 

In conclusion it should be emphasized however that even with all those theoretical un- 
certainties kept in mind, the problem of experimental observation of a new physical phe- 
nomenon, nuclear anapole moment, is a fascinating one. Moreover, if the theoretical results 
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of one-particle approximation, which are by themselves remarkably stable by nuclear stan- 
dards, will be supplemented by a serious treatment of many-body effects, those experimental 
investigations will give reliable quantitative information on P-odd nuclear forces. 
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Figure Captions 

Figure 1: First order correction 5R(r) to radial wave function of I/19/2 level. Dotted line 
corresponds to constant density and no spin-orbit potential (SR = rR(r)). Dash- 
dotted line refers to varying density and no spin-orbit potential. Dashed line 
refers to constant density spin-orbit potential included. Full line corresponds to 
full potential and varying density. 



Nucleus 133 CS 203,205 T1 209 Bi 207p b 

Analytical result || 
Its value at g p = 4 || 
Numerical result 
Its value at g p = 4 || 

1 

SI 

This work, analytical, g p = 4.5 
This work, numerical, g p = 4.5 

Table 1: Effective constants k for heavy nuclei 
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0.33 
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Nucleus 


\ r m )exp 


| r 2^2/3 


p\N + 3/2) 


41 Ca 20 


15.92(0.48) J22] 


10.27 


15.61 


45 Sc 21 
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16.11 


51 v 23 


12.89(0.72) H 


11.88 


16.79 


59 Co 27 


13.99(1.05) @ 


13.09 


17.62 



Table 2: Mean squared magnetic radii 
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Nucleus 133 CS 2 °3> 2 05 T1 209 Bi 207p b 



K, V 


0.317 


0.485 


0.376 


-0.099# n 


K S 2 ) 


0.301 


0.463 


0.349 


-0.095# n 


K. 3 > 


0.310 


0.497 


0.353 


-0.099# n 


«?„ 3) 


-0.088 


-0.120 


-0.126 





«L 3) 


0.014 


0.018 


0.019 


-0.006# n 


^4 1} 


0.003 


0.019 


-0.003 





54 2) 


0.01 


0.032 


0.014 





« 3) 


-0.019 


-0.055 


-0.019 





^ 3 ) 


0.048 


0.064 


0.070 





«J 3 ) 


-0.007 


-0.008 


-0.009 


0.004 


« *> 


0.268 


0.396 


0.309 


-0.105# n + 0.004 


« 2 ) 


0.258 


0.383 


0.298 


-O.lOl0 n + 0.004 


K 3 ) 


0.257 


0.396 


0.289 


-0.105# n + 0.004 



Table 3: Different contributions to the effective constant k. 
^ The potential parameters from Ref. |20| (see eq. fl69|)). 

2 ) The potential parameters from Ref. |25j (see eq. (|70|) ). 

3 ) The consistent LS-parameters (see eq. (|72|)). 
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